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Several exclusive B° decays are studied using a 1.86 fb -1 data sample collected at the T(5S) 
resonance with the Belle detector at the KEKB asymmetric energy e + e~ collider. In the B® — > 
D^tt + decay mode we find 10 B® candidates and measure the corresponding branching fraction. 
Combining the B° -> D s * ) ~n + , B° -» D^' p + , B° -> J/ip<t> and B° -> J/ipr) decay modes, a 
significant B° signal is observed. The ratio cr(e+e~ -> B*B*)/a(e + e~ -> B ( 3 * ) B S \* ) ) = (93±g ± 1)% 
is obtained at the T(5S) energy, indicating that B 3 meson production proceeds predominantly 
through the creation of B 3 B 3 pairs. The B 3 and B 3 meson masses are measured to be M(B 3 ) = 
(5370±l±3)MeV/c 2 andM(B s *) = (5418±1±3) MeV/c 2 . Upper limits on the B° -> 77, B 3 -> <jry, 
B 3 — > K + K~ and B° — > D s *^ + D s *^~ branching fractions are also reported. 

PACS numbers: 13.25.Gv, 13.25.Hw, 14.40.Gx, 14.40.Nd 



INTRODUCTION 

A considerable B° s production rate has been recently 
measured in e + e~ collisions at the energy of the T(5S) 
resonance [1, 2]. Thus, high luminosity e + e~ _B-factories 
have great potential for studies of exclusive B® decays. 
Although several B® decay channels have been recently 
observed by the Tevatron experiments [3, 4], a number 
of B° decay modes can be better measured at e+e~ col- 
liders running at the T(5S) energy. The detectors taking 
data at the T(5S) have many advantages in studies of 
B® decays, such as high photon and n° reconstruction 
efficiency trigger efficiency of almost 100% for hadronic 
modes and excellent charged kaon and pion identifica- 
tion. The possibility of partial reconstruction of specific 
decays and a model-independent determination of the 
number of initial B® mesons, which opens the possibility 
of precise absolute B° s branching fraction measurements, 
are additional advantages of B® studies at e + e~ colliders 
running at the T(5S). 

In this paper we report measurements of exclusive B® 
decays based on an T(5S) data sample of 1.86 fb _1 , col- 
lected with the Belle detector [5] at the KEKB asymmet- 
ric energy e+e~ collider [6]. This data sample is more 
than four times larger than the 0.42 fb _1 dataset col- 
lected at the T(5S) by the CLEO experiment in 2003 [7], 
where first evidence of exclusive B® decays at the T(5S) 
was found. 

We fully reconstruct six modes B® — > D~ir + , B® — ► 
D* s -tt+, B° s -> D~p+, B° s -> D* s - P +, B° - J/U and 
Bg — > J/iprj, which have large reconstruction efficiencies 
and are mediated by unsuppressed b — > c tree diagrams. 



Charge-conjugate modes are implicitly included every- 
where in this paper. To improve the statistical signif- 
icance of the Bg signal, these six modes are combined; 
the masses of the B® and B* mesons are determined from 
a common signal fit. 

In addition, we search for several rare B® decays: the 
penguin annihilation decay B® — ► 77, the electromag- 
netic b — > s penguin decay B° s — ► 07, and the hadronic 
b — > s penguin decay B° — > K + K~. Although the 
branching fractions for these decays are expected to be 
too small to be observed with this dataset, we can ob- 
tain useful upper limits. To date, only upper limits for 
the decays B® — > 77 [8] and B® — > (f>j [9] have been 
published. Within the Standard Model the B° — > 77 
decay is expected to proceed via a penguin annihilation 
diagram and to have a branching fraction in the range 
(0.5 — 1.0) x 10~ 6 [10, 11]. However, this decay is sen- 
sitive to some beyond-the-Standard Model (BSM) con- 
tributions and can be enhanced by one to two orders 
of magnitude in some BSM models [12, 13]. Although 
current measurements of the process B — ► X s j provide 
a more restrictive constraint for many BSM models, in 
these models the B® — > 77 process is more sensitive. 

The decay modes B® — > ^7 and B® — ► K + K~ are also 
mediated by penguin diagrams; these decays are natural 
processes in which to search for BSM physics [14-18]. 
The decay B® — ► K + K~ has been observed by CDF 
using a simultaneous multi-channel analysis [19], where 
overlapping signal peaks from the B® — ► K + K~ , B° — » 
K + n~ , B Q — » and B® — » K~n + decay modes 

were separated statistically in the fit. In this analysis the 
ratio (fj/fj) x B(B° S -+ K+R-)/B(B -+ K+^) = 
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0.46 ± 0.08 ± 0.07 was obtained, where (/J//J) is the 
ratio of production fractions of B° s and P° at Tevatron 
center-of-mass energy y/s = 1.96 TeV. 

We have also searched for the P" — > d{*^ + oi*' decay 
modes. These decay branching fractions are of special 
interest [20, 21]. These modes are expected to be pre- 
dominantly CP eigenstates and, because their branch- 
ing fractions are expected to be large, they should lead 
to a sizable lifetime difference between the CP-odd and 
CP-even P° mesons. Therefore within the SM frame- 
work the relative decay- width difference Ar^o /T b° can 

n S (*}+ S (*) — 

be obtained from measurement of the B" — > D s D s 
branching fractions. The first observation of the P° — > 
D+Dj decay has recently been published by the CDF 
collaboration [22]. 



BELLE DETECTOR AND EVENT SELECTION 

The Belle detector operates at KEKB [6] , an asymmet- 
ric energy double storage ring designed to collide 8 GeV 
electrons and 3.5 GeV positrons to produce T(4S) mesons 
with a boost of f3j — 0.425. In this analysis we use a 
data sample of 1.86 fb _1 taken at the T(5S) energy of 
^10869 MeV with the same boost. The experimental 
conditions for data taking at the T(5S) were identical to 
those for T(4S) or continuum running. 

The Belle detector is a general-purpose large-solid- 
angle magnetic spectrometer that consists of a silicon 
vertex detector, a central drift chamber (CDC), an array 
of aerogel threshold Cerenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) located inside a superconducting 
solenoidal coil with a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect 
K\ mesons and to identify muons (KLM). The detector 
is described in detail elsewhere [5] . A GEANT-based de- 
tailed simulation of the Belle detector is used to produce 
Monte Carlo event samples (MC) and determine efficien- 
cies. 

Charged tracks are required to have momenta greater 
than lOOMeV/c. Kaon and pion mass hypotheses are 
assigned based on a likelihood ratio £k/tt = £k/(£k + 
Ctt), obtained by combining information from the CDC 
(dE/dx), ACC, and TOF systems. We require £k/tt > 
0.6 (£k/tt < 0.6) for kaon (pion) candidates [23]. With 
these requirements, the identification efficiency for par- 
ticles used in this analysis varies from 86% to 91% (94% 
to 98%) for kaons (pions). A tighter kaon identification 
requirement Ck/tt > 0-8 is applied for the P° — > K + K~ 
decay, where the pion misidentification background is 
large. 

Electrons are identified combining information from 
the CDC (specific ionization dE/dx), the ACC, and the 



ECL (electromagnetic shower position, shape and en- 
ergy) [24]. Muons are identified by matching tracks to 
KLM hits and by using penetration depth information 
[25]. 

ECL clusters with a photon-like shape that are not 
associated with charged tracks are accepted as photon 
candidates. Primary candidate photons (7) that are 
used to reconstruct the P° — > (jry an d P° — » 77 de- 
cays are required to have proper bunch-crossing tim- 
ing and to lie within the acceptance of the ECL barrel 
(33° < 9 1 < 128°). To reduce the background from high- 
energy 7T° decays where the two daughter photons have 
merged into a single cluster in the calorimeter, the ECL 
energy deposition in a group of 3x3 cells is required to 
exceed 95% of that in the group of 5x5 cells around the 
maximum energy cell. The main background sources of 
high energy photons are 7r° — > 77 and 77—5-77 decays. 
To reduce these backgrounds, restrictions are imposed on 
the invariant mass of the candidate primary photon and 
any other photon (7') in the event. The primary pho- 
ton is rejected if 120MeV/c 2 < M(jj') < 145MeV/c 2 
and E 7 , > 30 MeV, or if 510McV/c 2 < M( 77 ') < 
570MeV/c 2 and P 7 / > 200 MeV. 

Neutral pion candidates are formed from pairs of pho- 
tons, each with energy greater than 150 MeV; the pho- 
tons must have an invariant mass within ±15MeV/c 2 
of the nominal ir° mass (i.e. <~ 3cr, where a ~ 5MeV/c 2 
is the 7T° mass resolution). A mass-constrained kine- 
matic fit is performed on the 7r° candidates to improve 
their energy resolution. We reconstruct r\ mesons only 
in the clean r\ — > 77 mode, requiring an invariant mass 
within ±20MeV/c 2 (<~ 2er) of the nominal 77 mass and 
photon energies larger than 50 MeV. Kg candidates are 
formed from 7r + 7r~ pairs with an invariant mass within 
±10MeV/c 2 (~3<r) of the nominal Kg mass and having 
a common vertex displaced from the interaction point by 
more than 0.1 cm in the plane perpendicular to the beam 
direction. 

The invariant mass for K* Q — ► K + n~ candidates 
is required to be within ±50MeV/c 2 of the nominal 
K*° mass; those of <f> — ► K + K~ candidates, within 
±12MeV/c 2 of the mass. A ±100MeV/c 2 mass win- 
dow is used to select p + — ► 7r + 7r° candidates. P7 mesons 
are reconstructed in the <fiir~ , K*°K~ and KgK~ de- 
cay channels; all candidates must have a mass within 
±12MeV/c 2 (~ 2.5er) of the nominal D~ mass. The 
P7 helicity angle distributions are expected to be pro- 
portional to cos 2 ^^ for pseudoscalar-vector final states; 
thus a |cos h *i\ > 0.25 requirement is applied for the 
D~ -> 4>n- and Dj -» K*°K~ decays. The helicity 
angle is defined as the angle between the directions 
of the K~ and Pj momenta in the <fi rest frame (or the 
directions of the ir~ and Pj momenta in the K*° rest 
frame in the case of K*°K~ decay). 

P*~ candidates are reconstructed in the P*~ — > P77 
mode; the measured D*~ and Dj mass difference is re- 
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quired to be within ±10MeV/c 2 of its nominal value. 
The invariant mass of candidate J/ip mesons is required 
to satisfy |M(//+ fi~) — mj/^\ < 30MeV/c 2 for the muon 
decay mode and satisfy — 100MeV/c 2 < M(e + e~) — 
m J/4) < 30MeV/c 2 for the electron decay mode, where 
mj/Tp is the nominal J/ip mass. 

B® decays are reconstructed in the following final 
states: D-jt+, D~p+, D*-^+, D*~ p+ , J/iJ><f>, J/ip V , 
Z?i*' + Z?i*'~, K + K~, (jry, and 77. The signals can 
be observed using two variables: the energy difference 
AE = Ega 4 — 2?bewn an d the beam-energy-constraincd 

mass M bc = ^/(E™J 2 - (p%f)^ where E™ and p™ 

are the energy and momentum of the B® candidate in the 
e + er ccntcr-of-mass (CM) system, and E^ m is the CM 
beam energy. The B® mesons can be produced in e + e~ 
collisions at the T(5S) energy via intermediate B*B* 7 
B*B° S , B° S B* S and B° S B° S channels, with B* s -» B° sl . 
These intermediate channels can be distinguished kine- 
matically in the M^ c and AE plane, where three well- 
separated Bg signal regions can be defined correspond- 
ing to the cases where both, only one, or neither of the 
B® mesons originate from a B* decay. The events ob- 
tained from MC simulation of the £?" — > D~ir + decay 
are shown in Fig. 1 for the intermediate T(5S) decay 
channels B*B*, B*B°, B°B* and B°B%. The signal re- 
gions are defined as ellipses corresponding to ±(2.0-2.5)<r 
(i.e. (95-98)% acceptance) resolution intervals in Mb c 
and AE. The signal events from the different interme- 
diate channels are well separated in the M\> c and AE 
plane. MC simulation shows that the separation between 
the channels in the M^ c projection is ~ 3 a or better 
for all studied B® decays. Elliptical regions do not de- 
scribe well the signal shape in the case of B® decays to 
the final states with photons or electrons, because the 
radiative energy losses result in a long tail on the left 
side of the signal AE distribution. In such decay modes 
the acceptance of the elliptical signal regions decreases 
to (70 — 80) %. A MC simulation indicates that the cor- 
relation between the M^ c and AE variables is small and 
can be neglected in this analysis. The numbers of events 
inside and outside these elliptical regions can be used to 
estimate the number of B® signal and background events. 

After all selections the dominant background is from 
e+e~ — > qq continuum events (q = u,d,s, or c). Topo- 
logically, B® events are expected to be spherical, whereas 
continuum events are expected to be jet-like. To suppress 
continuum background, we apply topological cuts. These 
were optimized using MC to model the signal and data 
outside the B® signal regions to estimate background. 
The ratio of the second to the zeroth Fox- Wolfram mo- 
ments [26] is required to be less than 0.3 for the high 
background 7r + , p + and K + K~ final states, 

less than 0.5 for the 77 final state (to increase the signal 
efficiency of such non-spherical B® decays) and less than 
0.4 for all the other final states. To suppress continuum 
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FIG. 1: The M bc and AE scatter plot for the B° DJn + 
decay obtained from the MC simulation. The ellipses show 
the signal regions for the intermediate B*B* (top elliptical 
region), B*B° and B°B* (middle elliptical region), and B°B° 
(bottom elliptical region) channels. 



further, the angle #* hr in the CM between the thrust axis 
of the particles forming the B® candidate and the thrust 
axis of all other particles in the event is used. We re- 
quire |cos0£ hr | < 0.9 for the low background final states 
with a J/ip, |cos#* hr | < 0.7 for the p + final states, 

|cos0* hr | < 0.6 for B° events reconstructed using the 
D~ — > K* Q K~ decay mode, |cos0* hr | < 0.5 for the very 
high background K + K~ final state, and |cos#* hr | < 0.8 
for all the other final states. 

More than one B° candidate per event can be selected. 
Using MC simulation we find that B° s decays to chan- 
nels with Dj or D* s ~ mesons can produce incorrect can- 
didates reconstructed in a cross-channel. Because the 
photon from the D*~ — > D~"f decay has a low energy, 
this photon can be removed from the B° reconstruction 
resulting in the replacement of an original D* s ~ by its 
daughter Dj or, conversely, an original Dj meson can 
be combined with a random photon to produce a false 
D*~ candidate. For example, every B° — > D*~n + de- 
cay will produce an incorrect B® — ► D~ir + candidate 
and ^37 % of B® — > D~ir + decays will produce incorrect 
B® — > -D*~7r+ candidates. Moreover, multiple candidates 
can be reconstructed in the B® — ► D*~tt + decay mode 
if the original photon from the D*~ decay is replaced 
by a random photon that satisfies the D*~ mass window 
requirement. The M\, c distribution of incorrectly recon- 
structed B° candidates has the same central value as the 
original signal, but the width is slightly larger. However, 
the AE distribution of these incorrectly reconstructed B® 
candidates is (200 — 300) MeV wide and shifted to neg- 
ative values if the correct photon is lost and to positive 
values when a random photon is added. 

We checked the effects of incorrectly reconstructed B® 
candidates on the results of the measurements reported 
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in this paper. Because of the large spread in the AE dis- 
tribution of the incorrectly reconstructed candidates and 
the low statistics used in this analysis, these effects are 
found to be small and are neglected; the corresponding 
uncertainties are included in the systematic error. We 
also checked other sources of multiple candidates in all 
studied decay modes and found that these effects can be 
neglected in the M^ c and AE measurements presented 
below. It should be noted that the MC efficiency calcu- 
lations also include multiple candidates and, therefore, a 
corresponding correction for this effect is applied. 



STUDY OF Bg -> D ( s * )_ tt 



D 



(*)- 



B° s 



J/ip<t> AND -»■ J/V>J7 DECAYS 



J/tp — ► e + e mode. As a cross-check, the branching 
fraction B{B° S -> J/ip<f>) = (0.9 ± 0.6 ± 0.2) x 10~ 3 is ob- 
tained for these two candidates, which agrees with the 
CDF measurement [29] within the large errors. The num- 
bers of B® candidates reconstructed in the D~ir + , D~ p + , 
£>*~7r + , D* s ~ p + , J/ip(j) an d J/ipV decay modes and lying 
in the signal region corresponding to the B*B* channel 
are listed in Table I. In addition, the numbers of events 
reconstructed in the three D~ decay modes are shown 
separately. 



TABLE I: The number of the B® candidates located within 
the elliptical signal region corresponding to the B*B* channel. 
The events reconstructed in the D~ — > <f>iT~ , D~ — > K*°K~ 
and — > K%K~ decay modes are listed separately. 



The M bc versus AE distribution for the B° -> D~-k + 
candidates is shown in Fig. 2a. Nine events are observed 
within the elliptical signal region corresponding to the 
B*B* pair production channel. Only one event is ob- 
served in the signal region for the B*B® + B°B* chan- 
nels, and no events are observed for the B° S B° S channel. 
Background outside the signal regions is small and corre- 
sponds to <~0.1 events in each of the three signal regions. 
The total number of bb events in the sample and the 
fraction of Bs B s events among all bb events at the 
T(5S) have been determined in [2] to be N% = (5.61 ± 
0.03 stat ± 0.29 syst ) x 10 5 and f s = (18.0 ±1.3 ± 3.2)%, 
respectively. We assume that 100% of B* mesons de- 
cay to the ground state B®. From the 10 observed 
events, the background estimate of 0.3 events, and the 
full reconstruction efficiency of (0.71 ± 0.10)% (inter- 
mediate branching fractions are included), we measure 
the branching fraction B(B° -> Djir+) = (0.68 ± 0.22 ± 
0.16)%. The systematic error includes the N^g and f s 
uncertainties and the uncertainty of ~ 14% in the re- 
construction efficiency, which is dominated by the uncer- 
tainty in the value of B(Dj — > </>7r~). This branching 
fraction is consistent with the value B(B® — » D~ir + ) = 
(0.38 ± 0.05 ± 0.14)% derived from a CDF measurement 
of B(B° S -» D-7r+)/B(B a -» D-TT+) [27] using the 2006 
PDG values of the B° — > D~ir + and D~ — > <j)ir~ branch- 
ing fractions [28]. 

Mb c and AE scatterplots are also obtained for the 

-,(*) 



Decay mode D s 



K*°K~ K%K~ Sum 



DV> p+ (Fig. 2c) 



B° s -» D*-n+ (Fig. 2b) and B° s 
decay modes. We observe four B° — > D* s ~-k + candidates 
and seven B° s — ► D s *^ p + candidates in the B*B* chan- 
nel, one B° s -> D { s * ) ~p+ candidate in the B*B° S + B° S B* 
channel, and no candidates in the B®B° S channel. 

The scatterplot in M bc and AE for the B° s — > J/ip(f) 
and B® — » J/xprj decays is shown in Fig. 2d. Two can- 
didates are reconstructed in the B° — > J/ip<fi mode and 
one candidate is reconstructed in the B° — > J/tpr] mode. 
One of the observed B® — > J/tp(f> candidates is recon- 
structed in the J/ip — > mode and one in the 



B° -> £>;- w + 

B° - D*~P + 
B° s -> 



Although the Mb c and AS signal resolutions are 
slightly different for different decay modes, for simplicity 
in Fig. 2 the same size elliptical signal regions are shown 
for all modes. Due to low statistics, a small variation of 
signal shape for different B° s decays can be neglected in 
the Mb c and AE measurements discussed below. It is, 
however, included in the systematic uncertainties. 

The six B® modes shown in Fig. 2 are combined 
to increase the statistical significance of the B® sig- 
nal. Distributions in AE arc obtained separately 
for events from three Mb c intervals, 5.408 GeV/c 2 < 

< 
< 



M bc < 5.429 GeV/c 2 (Fig. 3a), 5.384 GeV/c 2 < M bc 
5.405 GeV/c 2 (Fig. 3b) and 5.360 GeV/c 2 < M bc 
5.380 GeV/c 2 (Fig. 3c), corresponding to B® production 
proceeding through the B*B*, B* s B a s + B a s B* s or B°B° 
channels, respectively. 

Each of these three distributions is fitted with the sum 
of a Gaussian to describe the signal and a linear func- 
tion to describe the background. In the B*B* chan- 
nel (Fig. 3a), the width and the peak position are al- 
lowed to float, and their values <tae — (10.2 ± 1.9) MeV 
and (AE) = (-47.6 ± 2.6) MeV, respectively, are ob- 
tained from the fit. The width agrees with the value 
of ~ 12 MeV obtained from a MC simulation of the dom- 
inant B® — ► Dfn~ decay channel. Due to low statis- 
tics in the other two distributions, the peak positions 
and widths are fixed. The widths are taken from MC 
simulations. The peak position is fixed to zero for the 
B°B° channel and that for the B*B° S + B®B* channel is 
fixed to —23.8 MeV, which is half of the value obtained 
for the (AE) peak position in the B*B* channel. The 
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FIG. 2: The M bc and AE scatter plots for the B° -> D s ~7r+ (a), B° -> L>* _ 7T + (b), B° -> D^ )_ p+ (c) and B° -> and 
B° — > J/'i/'?? (d) decay modes. 
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FIG. 3: The AE distributions for B° 
5.408 GeV/c 2 < M bc < 5.429 GeV/c 2 , 
M bc < 5.405 GeV/c 2 and (c) 5.360 GeV/c 2 < M bc 
5.380 GeV/c 2 , corresponding to B° production through the 
B*B*, B*B® + B®B* and B°B° channels, respectively. 
Curves represent the results of the fits described in the text. 



fits yield 20.3 ± 4.8 events and 1.5 ± 2.0 events for the 
B*B* and B*B®+B®B* channels, respectively; no events 
are observed in the B®B® channel. From these numbers 
and approximately equal B® reconstruction efficiency in 
these three channels found in MC simulation, we ob- 



tain the ratio a(e + e 



B*B 



)/a(e+e 



Bi* ] Bi* ] ) = 



(93l£ ± 1)% at the T(5S) energy. The first uncertainty 
is statistical and the second uncertainty is systematic, 
dominated by uncertainties in the fit procedure. Poten- 
tial models predict the fraction of B*B* production to 
be around 70% [30-32]. 

The B* mass can be extracted from fit to the Mb c dis- 
tribution of the observed events in the B*B* channel. In 
this channel the Mh c variable, calculated from the for- 
mula Mbc = 



is equal, to a good 



approximation, to the mass of B* meson. This follows 
from the fact that the difference between the B° and B* 
momenta is statistically unbiased from zero and is smaller 
than the experimental resolution in B® momentum. Fig- 



ure 4 shows the Mb c distribution of the candidates in the 
range — 80MeV< AE < — 20MeV, where signal events 
from the B*B* S production channel are expected. We 
fit this distribution with the sum of a Gaussian to de- 
scribe the signal and a so-called ARGUS function [33] to 
describe the background. The fit yields a mass value of 
M(B* S ) = (5418±l±3)MeV/c 2 . The large systematic er- 
ror is dominated by the uncertainty in the collider beam 
energy calibration resulting in a e + e~ CM beam energy 
uncertainty of ~3MeV. The uncertainty of the method 
used to determine the M(B*) mass is estimated by MC 
simulation to be around 0.5MeV/c 2 . The uncertainty 
in the particle momenta measurements translated to the 
M(B*) mass uncertainty is also around 0.5MeV/c 2 . The 
observed width of the B* signal is (3.6±0.6) MeV/c 2 and 
agrees with the value obtained from the MC simulation, 
which assumes zero natural width and is dominated by 
the KEKB energy spread. The obtained B* mass is 1.8 a 
higher than the value measured recently by CLEO [34], 
M{B*) = (5411.7 ± 1.6±0.6)MeV/c 2 . 




(GeV/c*) 



FIG. 4: The B* mass distribution for events within the 
-80 MeV< AE < -20MeV interval, where the B° signal 
from the BgB* s channel is expected. Curve represents the 
result of the fit described in the text. 

Using the measured values of the B* mass and the 
energy difference (AE) = (—47.6 ±2.6) MeV, we can cal- 
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culate the B® mass. The value (AE) is the mean energy 
difference between the B* and B® mesons in the CM sys- 
tem and, in a good approximation, is equal to the mass 
difference of the B* and B® mesons. The photon energy 
in the B* — ► _B°7 decay is a constant in the B* rest frame, 
and the smearing due to the Lorentz transformation from 
the B* rest frame to the CM rest frame is small compared 
with the central value of the photon energy. Finally we 
obtain a mass value of M(B°) = (5370 ± 1 ± 3) MeV/c 2 . 
The second uncertainty in the B® mass value is the sys- 
tematic uncertainty dominated by the statistical uncer- 
tainty on the (AE) measurement, which will improve 
once more statistics become available. The uncertainty 
due to the collider beam energy calibration almost lineary 
affects both the M(B*) and (AE) values and nearly can- 
cels in the M(B Q S ) mass calculations. Other systematic 
uncertainties affecting the B® mass are similar to those in 
the B* mass measurement and are small. The obtained 
B® mass agrees well with the PDG value, M(B° S ) = 
(5369.6 ± 2.4) MeV/c 2 [28], and the most recent CDF 
measurement, M(B°) = (5366.01 ± 0.73 ± 0.33) MeV/c 2 
[35]. 



SEARCH FOR Bg -»■ 77, Bg ->■ <jyy, Bg -> K+K~, 
AND Bg -> D^+D^' DECAYS 

Distributions in M^ c and AE are also obtained for the 
reconstructed B° s — > 77 (Fig. 5a), B° s — » (jry (Fig. 5b), 
B° s -» K+K- (Fig. 5c) and B° -» d£* )+ d£* )_ (Fig. 5d) 
candidates. Only the £>° signal regions corresponding 
to the dominant B*B* channel are considered for the 
searches reported here. These regions are wider for the 
B® — > 07 and £>° — > 77 decays, where energy losses due 
to photon radiation lead to a large tail at lower values of 
AE. The signal region for the B® — > Dg*'" 1 "^*' - search 
is slightly smaller, because of the kinematics of the decay 
to two heavy particles. The shapes of the signal regions 
for these decays are optimized from the MC simulation. 

To avoid multiple D* and D s cross-channel candi- 
dates, only one candidate per event is selected in the 
Bg — > Dl*^ + Z)g*^ analysis, where the number of multi- 
ple candidates can be rather large. The candidate with 
the Mbc value closest to the nominal M(B*) value mea- 
sured above is chosen. No significant signals are observed 
in any of the distributions shown in Fig. 5. However 
one Bg — ► 07 event, two B® — > K + K~ events and 
one Bg — ► D* + D~ event lie within the signal regions, 
whereas backgrounds outside the signal regions are not 
large. The numbers of events within the signal regions, 
the estimated background contributions, the efficiencies, 
and the upper limits for the corresponding B® branching 
fractions are listed in Table II. For comparison, the pre- 
viously published upper limits and branching fractions 
are also shown. The numbers of events and the upper 



limits are obtained using only the B s signal region cor- 
responding to the B*B* channel. The upper limits are 
obtained using the Feldman-Cousins method [36], and a 
small correction due to systematic uncertainties is ap- 
plied. The efficiencies are determined from the MC sim- 
ulation. The number of initial B* S B* S pairs is obtained 
by multiplying the number of Bs B s pairs measured 
in the inclusive analysis [2] by the production ratio of 
B*B* pairs to all B^ B^ pairs obtained in this analy- 
sis. We calculated the previous B° — > DfDj branch- 
ing fraction listed in Table II using the measurement 
B(B° -» D+D-)/B(B° -» D+D~) = 1.67 ± 0.41 ± 0.47 
from CDF [22] and the B(B° D+D-) value from the 
PDG [28]. 

The upper limit obtained for the decay B® — > 77 is 
about three times smaller than the most restrictive pub- 
lished limit [8]. However, it is still two orders of mag- 
nitude above SM predictions [10, 11]. The upper limit 
obtained for B° s — ► <jry is about a factor ten larger than 
the theoretically expected branching fraction [37]. The 
upper limit obtained for the B® — > K + K~ decay is an or- 
der of magnitude larger than the value measured by CDF 
[19]. For SM branching fractions, statistically significant 
signals of ~10 events can be obtained for the B° — > 4>-f 
and B® — > K + K~ modes in a ~30fb _1 dataset on the 
T(5S). 

The upper limits obtained for B° s -> D s * )+ D s * } de- 
cays are of special interest because the D^*^ states 
are expected to be dominantly CP eigenstates. Assum- 
ing that the branching fractions for the DfDj, D* + D~ , 
D+D* s ~ and D* + D* s ~ final states are each in the range 
(1-3)%, we expect about 5-10 events in each of these 
four channels with statistics of ~30fb _1 . Within the SM 
framework such measurements can provide an important 
constraint on the value of Ar B o/r B o [20, 21]. 

CONCLUSIONS 

Several exclusive B° decays are reconstructed using 
1.86 fb -1 of data taken at the T(5S) resonance with the 
Belle detector at the KEKB asymmetric energy e + e~ col- 
lider. 

R° signals are found in six decay modes: B° s — > 
Z?i* K 7r+, B° s -» Di* ) ~p+, B° s -» J/ij<t> and B° s -» 
J/ipff. The branching fraction B(B® — > D~ir + ) = 
(0.68 ± 0.22 ± 0.16)% is measured. Combining the stud- 
ied six channels, we observe a significant B° s signal and 
obtain the masses M(B° S ) = (5370 ± 1 ± 3)MeV/c 2 
and M(B*) = (5418_± 1 ± 3)MeV/c 2 . B° produc- 
tion through the B*B* channel is found to dominate 
over other B^B^ channels; the ratio er(e + e~ — » 
B*B*)/a(e+e- -» B^B^) = (93^±1)% is measured. 
These results are in agreement with CLEO measurements 
[7]- 



8 



5.5 



U 

> 5.4 
CD 

o 



5.3 




5.2 



-0.4 -0.2 



i 

0. 



(b) 



0.2 



-0.4 -0.2 



. i i 
0. 



' l 1 1 1 l 1 1 1 l 

(c) 



0.2 



-0.1 



0. 



AE (GeV) 



(d) 



0.1 



-0.2 -0.1 



i 

0. 



0.1 



0.2 



FIG. 5: The scatter plots in M hc and AE for the B° 77 (a), B° -> <t>j (b), B° -»■ K+K~ (c) and B° -> D^D^ (d) 
decay modes. In the latter case, the signal event is reconstructed in the B® — > D* + D~ decay mode, while the three background 
events are reconstructed in the B° — ► Dj D~ decay mode. The ellipses indicate the B® signal regions for the B*B* channel. 



TABLE II: The number of events in the signal region (Yield), the estimated background contribution (Bkg.), the efficiencies 
(Eff.), the 90% C.L. upper limits derived in this analysis (Belle upper limit) and previously published upper limits or branching 
fractions (Previous UL/BF) for the B° -> 77, B° -> 07, B° -> K + K~ and B° 3 -> D ( s * )+ D^' decay modes. 



Decay mode 


Yield 


Bkg. 


Eff. 


Belle 


Previous 




(events) 


(events) 


(%) 


upper limit 


UL/BF 


B° s -> 77 





0.5 


20.0 


<0.53 xlO" 4 


<1.48xl0" 4 [8] 


B s ° -> 07 


1 


0.15 


5.9 


<3.9 xl0~ 4 


<1.2 xl0~ 4 [9] 


Bg -f i^A^ 


2 


0.16 


9.8 


<3.1 xl0~ 4 


(3.30 ±0.57 ±0.67) x 10" 5 [19] 


B s ° - D+D7 





0.02 


0.020 


<6.7% 


(1.09 ± 0.27 ± 0.47) % [22] 


B s ° - D*+D~ 


1 


0.01 


0.0099 


< 12.1 % 




B ° -> D*+D*~ 





<0.01 


0.0052 


<25.7% 





We have also searched for — ► 77, B° s — > <j>j, 
B° s -» K+K- and B° s -» Di* )+ Z?i* K decay modes and 
set upper limits on their branching fractions. The upper 
limit on B® — > 77 is three times more restrictive than 
the best existing limit. The background levels in these 
decays are low, indicating that the sensitivity of future 
studies of these decays with larger statistics will not be 
limited by backgrounds. We expect that significant sig- 
nals for B° s -> K+K~, B° s -> 7 and B° s -» Di* )+ Di* } ~ 
decays can be observed in ^30fb _1 of data. With such 
statistics the upper limit for the B® — > 77 decay should 
provide an important constraint on some BSM models. 
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